Quantum Hydrodynamics, Quantum Benjamin-Ono Equation, and Calogero Model 
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Collective field theory for Calogero model represents particles with fractional statistics in terms 
of hydrodynamic modes - density and velocity fields. We show that the quantum hydrodynamics of 
this model can be written as a single evolution equation on a real holomorphic Bose field - quantum 
integrable Benjamin-Ono equation, ft renders tools of integrable systems to studies of nonlinear 
dynamics of ID quantum liquids. 



1. Introduction. Calogero-Sutherland model occupies 
a special place in ID quantum physics. A singular in- 
teraction (an inverse square potential) represents a frac- 
tional statistics of particles (for a review see e.g., 0]). 
The eigenfunctions of the model are neither symmet- 
ric nor antisymmetric. They have a form (A(x)) x J(x), 
where A is a Vandermonde determinant and J are sym- 
metric polynomials (Jack polynomials). Although the 
problem does not possess a conformal invariance, it is 
intrinsically related to a boundary conformal field the- 
ory (CFT) 0] with a central charge c — 1 — 24a 2 ,, where 
2ao = \/\ — 1/ VA. Here we focus on a rational degener- 
ation of the model - the Calogero model (CM) 
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One might also add a harmonic potential \lo 2 ^ x\ to 
prevent particles from running to infinity. Most of the 
formulas below are simple to generalize to the chiral sec- 
tor of the trigonometric version of the model. 

A productive approach to Calogero model in the limit 
of an infinite number of particles is a collective field the- 
ory developed by Sakita and Jevicki in Refs. 0, Q and 
extended to CM in Refs. [HQ In this approach an 
entire quantum theory with all the richness of an opera- 
tor content, is represented as a quantum dissipationless 
hydrodynamics j6jj, i.e., solely in terms of the density 
p(x) — J2i 3( x ~ x i) an d the velocity operators 

\p(x),v(y)] = -iS'(x - y). 

Hydrodynamic approach is especially useful if one is in- 
terested in the evolution of smooth wave packets, larger 
than an inter-particle distance. Original particles apear 
in this approach as solitons of nonlinear fields. 

In this letter we show that the quantum hydrody- 
namics (aka bosonization, or collective field theory) of 
Calogero model is equivalent to a quantum version of an 
integrable Benjamin-Ono equation on the double (QBO). 
The latter is a minor generalization of the conventional 
Benjamin-Ono equation Q30) arising in the hydrody- 
namics of stratified fluids 7] . A connection of the tradi- 
tional BO to a complexified version of a classical Calogero 



model is known |8|. A pair of coupled classical BO-type 
equations has been obtained by A. Jevicki for a long-wave 
description of free fermions (A = 1) p|. 

2. Quantum hydrodynamics and non-linear bosoniza- 
tion. Collective or hydrodynamic approach is a nonlin- 
ear version of the bosonization procedure [9j - a popular- 
method in low dimensional physics. In this approach, a 
fermionic spectrum is not linearized at the Fermi level. 
Therefore, an asymmetry of the particle-hole spectrum - 
a dispersion of hydrodynamics modes - is not neglected. 
The price is a nonlinear and dispersive character of the 
hydrodynamics. The gain is a realm of nonlinear phe- 
nomena missed in a linear approximation. Before pro- 
ceeding to the Quantum BO equation on the double 
(QBO), we start with the quantum hydrodynamics of 
free fermions. 

Quantum hydrodynamics of free fermions with 
quadratic dispersion (A = 1) can be obtained using con- 
ventional "bosonization". Writing a fermion tp( x ) — 



as a superposition of left and right chiral 



fields [d x ifL,R{x), </?i,it(0)] = =pi5(x) one obtains 



2ir(n + l)^(-id x ) n i> ~: (d x <p R ) n+1 - (d x cp L ) 



n+l 



(2) 



modulo full derivatives. Here and below we assume nor- 
mal ordering. The first three moments are density p = 
^d x (ip R -(p L ), current j = pv = ■^((d x (p R ) 2 -(d x (p L ) 2 ), 
and Hamiltonian density 



H 
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{{d xWR f - (d xVL f) = \pv 2 + yp 3 . (3) 



The operator v — \d x ((p R + tpi,) is a velocity opera- 
tor. The meaning of the terms in the energy J2J is 
clear. The "kinetic energy" reflects the Galilean invari- 
ance of the "potential energy" reflects the fermionic 
statistics: (7r 2 /6)p 3 is an energy density of a Fermi 
gas. This Hamiltonian takes into account all long-wave 
correlation properties of fermions, missing the physics 
at 2k R (e.g., Friedel oscillations), where the chiral sec- 
tors interact. Linearization at the equilibrium density 
po = kp/ir leads to a familiar linear bosonization H w 
f {{d xVR ) 2 + [d xVL ) 2 ) = ^(v 2 + 7T 2 (p-p ) 2 ). This ap- 
proximation misses an interaction and dispersion of hy- 
drodynamics modes, essential in many applications. 
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The Hamiltonian (J3J) generates operator equations of 
quantum hydrodynamics 6] 



continuity equation 
Euler equation 



p + d x (pv) = 0, (4) 
v + d x {\v 2 + w) = 0, (5) 



where w — 5(pe)/6p is the quantum enthalpy and e = 
(ir 2 /6)p 2 is the energy of a filled Fermi sea per particle. 

3. Holomorphic Bose field. Both continuity equa- 
tion J3J and Euler equation JSJ) for free fermions can be 
compactly written in terms of a single holomorphic field 
u(z,t). It is a real field obeying the Schwarz reflection 
symmetry u(z) — u(z). This field glues together left 
and right components of the Bose field. Let us pass to 




FIG. 1: The double. Left and right sectors are glued together 
along the branching cut (equator) - a segment of the real axis 
supporting a particle density. Excitations are symmetrically 
located in a complex plane. 



an imaginary time t — > —it, v — > iv, and treat ID co- 
ordinate real axis of a complex plane z. Then 
right and left currents ±d x (pn.L — irp ± iv are com- 
plex conjugates and may be treated as boundary val- 
ues of an analytic field in a domain adjacent to the real 
axis u(x ±i0,t) = v ± inp = —id x <pB..L- Then the hy- 
drodynamic equations I|4I5[) are seen as boundary values 
(z = x ± i 0) of the quantum Hopf-Burgers equation 



dtu + ud z u = 0. 



(G) 



This is a nonlinear bosonization of free fermions. 

The analytical structure allowing to treat the right cur- 
rent as an analytical continuation of the left one across 
the boundary (the real axis x) extends to the Calogero 
model. There, as we argue below, the boundary values 
of the field u are 

u{x ± i 0) = (l/\/X)(v ± iXirp) + a d x logp, (7) 

The field u consists of positive and negative parts 

u = u~ (z) + u + (z). 

The negative part is a Cauchy transform of the density 

p{x) 



(z) = ^Jdx- 



It is an analytic function in a plane cut by a seg- 
ment of a real axis supporting the density. Its bound- 
ary values u~(x ± iO) = ii\f\{p H ± ip), where p H = 
(l/7r)P J dx' p(x')/(x' — x) is a Hilbert transform of the 
density. A positive part - u + (z) is an analytic continu- 
ation of the real function u + (x) = (l/yX(v — Xnp H ) + 
aod x \ogp from the real axis. It is analytic only in the 
vicinity of the real axis, having singularities elsewhere. 

Let us mark two points, one inside, another outside 
of the domain, where both u + and u~ are analytic, and 
expand the fields around these points. For simplicity let 
us talk about the rational case. The natural points are 
z = and z = oo. We have 



n>0 



n>l 
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The modes a n are real and Hermitian. Positive modes 
are moments of density (collective variables) a n — 
\f\ J x n p(x)dx and negative modes are canonically con- 
jugated to positive ones [a„, a_ m ] = —nb n ^ m . The zero- 
mode is the number of particles ao = yXN. It docs 
not change in time. The field u{z) is a current of a real 
canonic Bose field ip = ip + + <p~ , where d z ip ± (z) — ^(z). 

Let us notice a difference between the Bose field (JHJ 
and a Bose field used in conventional bosonization of free 
fermions with a linearized spectrum. In the former, the 
time t is not a part of a z-plane, but describes an evo- 
lution of a real holomorphic field u{z,t). In the latter, 
time and space form a complex coordinate x + it of the 
radial quantization. Complex z-plane has an analytical 
structure similar to ones used in boundary CFT and is 
sometime referred to as (Schottky) double. 

4- Quantum Benjamin- Ono equation is defined as 



(), ,t - i), ( iu 2 + a d z (u + - u ) 



0. 



(9) 



where we rescaled time t — > tj\f\. As the Hopf-Burgers 
equation JfjJ, it is a quantum equation acting in the Fock 
space spanned by modes a n . QBO is integrable. It has 
infinite number of conserved integrals. The first four are 
u itself, momentum density T(z) = u 2 /2+iaod(u + — u~), 
the Hamiltonian 



H = 



dz 
Aiti 



u + actud z (u + — u ) 



(10) 



and the fourth integral 
1 3a 



h = 



dz 
2ttI 



u 2 d z (u+ - u~) + 2al(d z uf 



(11) 

Here the integration goes around the real axis. Below we 
show that HOIlOfl are equivalent to the CM (JIJ . 

The classical BO equation originally appeared in hy- 
drodynamics of stratified fluids. There the real harmonic 
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function u{x,y) on a half plane evolves followinga time- 
dependent boundary condition on the real axis [7J 



dtu + ud x ii + d„u 



0. 



fl2) 



The Hilbert transform u H = (l/ir)P J dx' u(x')/{x' — x) 
can be seen as u H — i(u + — u~) with respect to 
the one marked point at infinity. Contrary, our BO- 
equation has two the marked points z = and z = 
oo. The term (u + — u~) = —iuq in (0 can be 
seen as a Hilbert transform with respect to the con- 
tour C embracing a support of the density (real axis) 
S C) — —P § c u ^ z _l dz' . This difference is not an 
obstacle to an almost automatic extension of known so- 
lutions of the traditional BO equation. Remarkably, this 
extension carries over to the quantum case as well. 

5. Quantum Hirota equation. Typically, integrable 
equations can be written in the bilinear (Hirota's) form 
(see for a review [lOL 1 1 1| ) . Apparently, bilinear forms ex- 
ist also for quantum integrable equations. Let us intro- 
duce a vertex operator ty a — exp(— cr<p) : where a = ^f~\ 
or <7 = — l/VA, and consider its positive and negative 
parts 



-aip + 



The bilinear quantum Hirota equation 



D t i/j- • V> H 



(13) 



(14) 



is an equivalent form of QBO ©■ Hirota's derivatives 
are defined as D x a ■ b — (d x — d x ') n a(x)b(x')\ x — x >. The 
proof is straightforward. 

The Hirota equation (|14f> is a compact form of OPE 
(J2J, extended for Calogero model. The resemblance to 
the Schrodinger equation further clarifies the origin of 
QBO. It is especially transparent for free fermions A = 1. 
There, the vertex operator 'f a is a bosonized form of a 
fermion creation [a = 1) and anihilation (<r = —1) oper- 
ators. The Schrodinger equation for free fermions = 
\d 2 ^ transforms to Hirota equation 114(1 by an elemen- 
tary algebra. A known duality between particles and 
holes a n — > a_„ under a transformation \f\ — > — 1/vA is 
especially transparent in the bilinear form. 

6. Quantum hydrodynamics and QBO-equation. Re- 
lation to hydrodynamics is obtained by specializing BO- 
equation on the real axis. Computing the jump (discon- 
tinuity) of the l.h.s. of © across the real axis we obtain 
the continuity equation I0J with the current j = pv = 
(T(x + iO)-T(x-iO))/{2mVX). It gives the announced 
connection J7|) between the field u and velocity. The sum 
of BO equations taken on two sides of the real axis gives 
the Euler equation J^l with the energy density 

e = 1(tt\ P ) 2 + ^(d x lnp) 2 +na \d xPH . (15) 
o 2 

Eqs. (|4I5I15|) have been derived in Ref. 5] within the 
collective field theory approach. 



7. Pole dynamics of quantum Benjamin- Ono equa- 
tion and Calogero model. Let us now show that QBO is 
equivalent to the Calogero model. 

An important class of solutions of classical integrable 
nonlinear equations is given by a moving pole ansatz. 
These solutions have a form u = X)j=i f( z ^ x i(t))i where 
f(x) is an elliptic function, or its trigonometric or ratio- 
nal degeneration (see, e.g., (tJ). Then nonlinear equation 
is reduced to a certain, also integrable, many body prob- 
lem of N particles with coordinates Xi. In most cases 
it is the Calogero model and its decendents. In Ref. 
H similar fact has been noticed for the classical BO- 
equation 1(12(1 • The complex poles Xi (t) of a pole ansatz 
u(x,t) = - J2 l ( x 2 = x - - ^31") move according to a com- 
plexified classical Calogero equation 



2£ 

3,3^ 



X 2 



(Xi Xj) 



(16) 



under the condition that Imx, > at an initial time. 

Our version ij2*T|) of the classical BO equation produces 
a genuine Calogero model where coordinates of particles 
(poles) are real. Remarkably, the pole ansatz carries over 
to the quantum case. Pole ansatz is perhaps a simplest 
way to see the equivalence between CM and QBO. 

The algebra is straightforward and can be carried out 
for both QBO and quantum Hirota equation. We briefly 
describe the main steps starting from the QBO ©■ 

We look for a solution in the form 



JV 



i 



J X Z- Xi(t) ' 



= Real 



(17) 



and show that the dynamics of the poles xi obey quan- 
tum Calogero model with N particles. The residue in 
()17(l is controlled by the asymptote at infinity u~(z) ~ 
(—>/XN)/z. At the pole ansatz the vertex operator 1(13(1 
at a = — l/\/A has simple zeros at positions of particles 

= Uj( z - x 3)- 

We plug ((17(1 into the negative part of @ and assume 
[u + (xi),Xj] — —Sij/VX, which yield canonical commu- 
tation relations for u as N — > oo. Using formulas 

(u-u+)+ = - ]T\ T^-( u+ ( z ) - u+ ( x i)) and («~« + )~ = 

E. — u + (xa ). we obtain that a rational function 



k,k^j 
1 



Z tJU'H Z X n 



y/\(z-Xj) 3 (Z-Xj)'- 



(Xj - x k ){z - Xj ) 2 
u+( Xj ) = (18) 



vanishes. Both second and third order poles vanish if 
Xj = pj with [p*, afj] = —Sij/y/X and 



= u+ ( x 3) - ^ E 



1 



k,k=£j 



Xj Xk 



(19) 
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Eq. (|19fl is a finite N version of the equation J7J • Differ- 
entiating (|19|l over time and using the positive part of 10 , 
we obtain after some work Ea. (|16|l with A 2 — > A(A — 1) 
(where time is scaled back to the original real time). 

Few comments are in order: (i) The pole ansatz l|17f) 
represents positive modes a„ = \/A J2j x ^ ■ Negative 
modes are represented in terms of momenta and coordi- 
nates of particles through l|19(l . Canonical commutation 
relations between positive and negative modes appear 
only in the field theory limit. For a finite number of par- 
ticles they are replaced by the conditions [u + (xi), Xj] = 
—Sij/VX used in the cancelation of poles; (ii) Simi- 
lar manipulations in the classical case produce Calogero 
model with interaction X 2 (xi — x.j)~ 2 . A quantum shift 
A 2 — ► A(A — 1) comes from pulling momentum operator 
Pj to the right in the first term of (|18|l : (iii) A trigono- 
metric ansatz u~(z, t) = —yXk cot(fc(z — Xi(t)) leads 
to the Calogero-Sutherland (trigonometric) model. 

8. Solitons and integrability Let us now consider the 
large N limit and assume that there is a non-zero back- 
ground density of particles. Then, alternatively, one 
looks for the pole ansatz for u + with a finite number 
of poles symmetrically located with respect to the real 
axis (see FigQ 

u^,t) = Vx±( 7 ^ i + ^). (20) 

In this case the vertex operator (|13[) at a = — 1/vA is 
^jj~{z) = Yl^z — z^ and has zeros at positions of poles. 
|2| These poles also move according to quantum Calogero 
model. They correspond to excitations of a system with a 
continuous density. Accordingly u~{z) has a branch cut 
on the real axis. Solutions (states) with well-separated 
poles Zi of u + (z, t) are quantum analogs of solitons known 
in classical nonlinear equations. A semiclassical one- 
soliton solution of the Euler equation JSJ) studied in Ref. 
0, 0, El , corresponds to a single pole in (|20|l and can 
be found from Q). 

Existence of solitons signals that the QBO is inte- 
grable. Indeed, the known Lax representation for the 
classical BO equation lfl2^l extends to the quantum BO 
equation on the double (0 with minor changes. In partic- 
ular, the pole ansatz converts Lax operators of the quan- 
tum BO equation to known Lax operators of quantum 
Calogero model. The same is true for the conserved inte- 
grals, Baklund transformation, etc. In a standard way 
these tools provide multi-soliton, periodic, and shock- 
wave solutions for quantum hydrodynamics. 

10. Semiclassical limit of BO equation and nonlinear 
transport in quantum systems. Classical BO equation 
may be useful in studying nonlinear effects in transport 
in low-dimensional quantum systems. It describes a semi- 
classical evolution of long wave packets of density. The 
semiclassical limit is subtle and requires certain care. In 
nonlinear quantum equation @ a simple replacement of 



a quantum field it, by classical fields is incorrect. The 
correct result is obtained by shifting ao — > V~\/2 in © 

d t u + d z Qm 2 + (V\/2)d z (u + - u'^j = 0. (21) 

In particular, the dispersive (last) term in (|21|l does not 
vanish for free fermions at A = 1, while it vanishes in 
the quantum version (jfJJl. The equation l|21|) (or rather 
a corresponding classical action) must be understood as 
a leading term of the gradient expansion of the effec- 
tive action, with the first order correction due to zero 
fluctuations. It can be summurized as (u 2 ) — (u) 2 
V A d z (u + — u~ ) , where we average over the coherent state 
with a given density. The correction ammounts to the 
shift A(A — 1) — * A 2 . One can also obtain this quantum 
correction by comparing the classical i|16[) and quantum 
Q CM obtained as pole ansatzes of classical and quan- 
tum BO equations. The classical pole asatz Ijl7[l leads 
to the classical Calogero equation Qlfijl only when it is 
plugged into H21(l rather than into a naive classical ver- 
sion of (0 • The classical dynamics of long wave packets 
depends on A only through scales. It is described by 
the classical BO on the double, and is the same for free 
fermions and "anyons" . 

The appearance of the dispersion term (the last term 
in JHJ) for free fermions has been pointed out by A. 
Jevicki 4]. It has far reaching consequences. As it is 
well known, a dispersionless limit of non-linear equations 
is singular |l4j |. Almost any smooth initial configura- 
tion u(x, t) evolving according to a classical Hopf-Burgers 
equation © (a dispersionless limit of J2D) develops an 
unphysical shock wave in finite time. The dispersion be- 
ing initially small becomes a crucial factor and can not 
be neglected. A discussion of the role of dispersion in 
electronic transport, and a number of other interesting 
topics, like quantum integrability, relation to perturbed 
conformal field theory etc., is outside of the scope of this 
letter. We plan to address some of these issues elsewhere. 
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